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The existence of a Goldstone boson is demonstrated in 
chiral random matrix theory. After determining the effective 
coupling and calculating the scalar and pseudoscalar propaga- 
tors, a random phase approximation summation reveals the 
massless pion and massive sigma modes expected whenever 
chiral symmetry is spontaneously broken. 

PACS: 11.30.Rd, 11.30.Qc, 12.38.Lg. 

Increasing evidence supports the applicability of chi- 
ral random matrix theory (%RMT) fij for describing the 
spontaneous breaking of chiral symmetry in quantum 
chromodynamics (QCD). (For a review, see Ref. 0.) 
Using the Banks-Casher relation ||, the vacuum chi- 
ral condensate, S, can be related to the spectral density 
of the Dirac operator in the vicinity of zero eigenvalue. 
This model has been extended to describe chiral symme- 
try restoration, for both finite temperature and baryon 
chemical potential Q], in what amounts to a mean field 
treatment. Although these efforts have led to a schematic 
phase diagram for QCD as a function of T, /i, and the 
current quark mass, m, other familiar vacuum phenom- 
ena associated with chiral symmetry breaking remain to 
be obtained. In this letter, we consider quark-antiquark 
correlation functions in order to recover the Goldstone 
mode. This is a ubiquitous feature associated with the 
spontaneous breaking of a continuous global symmetry 
which is expected whenever the symmetry group is of 
order greater than one. In xRMT, SV L (N f )xS\J R (N f ) 
chiral symmetry is broken into SUy(iVf) as is believed 
to occur in simplified QCD with two massless flavors. 
Along with the chiral condensate, (■tpip)o, we should ex- 
pect xRMT to provide us with massless pseudoscalar 
(i.e., the pion) and massive scalar (i.e., the sigma) ex- 
citations. 

In order to extract this feature from a model of chiral 
symmetry breaking at the quark level, a resummation in 
the spirit of the random phase approximation (RPA) can 
be performed to construct the current-current correlation 
function in a given channel. This is shown schematically 
in Fig. 1. Evidently, this procedure requires two model- 
dependent quantities: an effective interaction between 
quarks and a correlator describing the propagation of a 
quark-antiquark pair in the gluonic background. If we 
denote the coupling as g and the "bare" qq propagator 
as IIo, the full or "dressed" correlation function, II, can 
be written as an infinite sum. Each term in this sum 
is an alternating string of coupling constants and bare 
propagators as indicated in Fig. 1. This geometric series 



is readily written as: 
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In principle, both g and IIo can depend on the momentum 
carried by the quark-antiquark pair in a model-dependent 
manner. We now turn to the particular case of xRMT. 




FIG. 1. The qq propagator in the random phase approxi- 
mation. 

To determine II, we assume a factorization of the ele- 
ments g and IIo , in which case its solution takes the RPA- 
like form of Eq. ([!]) and Fig. 1. The coupling constant is 
computed as the average strength of the four-quark inter- 
action induced by the random matrix background. Next, 
IIo is schematically associated with the two-eigenvalue 
correlation function, calculable in %RMT, and evaluated 
in the appropriate limits. This is the zero-dimensional 
analogue of a loop integral, in that it encodes the effects 
of the chiral condensate on a quark-antiquark pair. 

The effective coupling, g, is easily determined by aver- 
aging over the random matrix elements in the partition 
function of xRMT in order to generate an effective quark 
interaction. This partition function can be written as 
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where the W are nxn complex matrices and where tp and 
tjy are row and column vectors with 2n components rep- 
resenting two-component chiral spinors. Implicit in this 
form is the choice of a Gaussian unitary ensemble (GUE), 
appropriate for the SU(3) gauge group, with variance E. 
The choice of a Gaussian potential is not mandatory in 
chiral random matrix theory since the microscopic spec- 
tral correlators have been shown to be universal for a 
wide class of potentials H||. 

The integral over W in Eq. (g) can be performed to 
obtain 
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from which one immediately finds an induced four-quark 
interaction. After normalizing for N, the effective cou- 
pling strength is g — (2JV/E 2 ) -1 . This coupling is con- 
stant, as is to be expected from the mean-field nature of 
random matrix theory. 

The second ingredient, the bare current-current corre- 
lation function, has the form 
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Tr [iS(p)TiS(p - q)T] 
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Here, S(q) is the fermion propagator, T is one of the 16 
Dirac basis matrices which generate the Clifford algebra, 
and q is the momentum carried by the quark-antiquark 
pair. The two relevant correlators here are the identity 
matrix T = 1 and T = 75 appropriate for the scalar 
and pseudoscalar channels, respectively. We will be con- 
cerned only with results valid in the q — > limit of zero 
momentum transfer. 

Fortunately, the bare propagators of Eq. (^) are di- 
rectly related to the function, p(Xi,X2), describing the 
correlation between two eigenvalues of the Dirac opera- 
tor. This spectral correlator is accessible in xRMT and, 
moreover, known to be universal in the microscopic limit. 
The calculation of spectral correlators is usually techni- 
cally involved, but it simplifies significantly in random 
matrix theory through the use of orthogonal polynomial 
techniques . The desired spectral correlator is 
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for eigenvalues Ai and A2 generated by averaging over 
N x N matrices H with elements drawn from a random 
Gaussian distribution. For problems with chiral symme- 
try, H has the form 
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The scalar and pseudoscalar correlation functions of 
Eq. (Ji|) can be constructed from this spectral correla- 
tor. We will be interested in taking limits N — > 00 and 
A = (Ai — A2) — > simultaneously. This requires first 
fixing A ~ 1/N and then taking TV — ► 00. Since the 
effective coupling of Eq. (|^) does not depend on the mo- 
mentum transferred, we are interested ultimately in the 
limit in which A — > 0. The limits cannot be interchanged, 
for their ordering is crucial in obtaining the physically 
relevant quantities from xRMT B. 



To extract the A dependence in the described limit we 
use the microscopic spectral correlation function, 



p M (x,y)= lira p( x/NY, , y/NT, ) . 



(6) 



The microscopic limit is taken in order to retain infor- 
mation on the scale of the spacing between neighboring 
eigenvalues, the scale of 6 = ANY which corresponds 
to the momentum q. At the same time the eigenval- 
ues themselves are rescaled into microscopic variables 
x = XNY and the limit A — > is replaced by x ^> 1. 
Thus, we are considering two eigenvalues which are close 
on the macroscopic scale, but have a difference much 
greater than the average spacing between microscopic, 
unfolded eigenvalues. 

In these limits the real parts of the traces of Eq. (|5|) 
vanish, being odd in A, and we need only consider the 
contributions from the imaginary parts. These may be 
found by taking the discontinuities across the real axes 
in the limit of zero current masses to, to'. For a single 
fermion species, 

lim Il(g 2 ) = lim lim Disc| m= jA,m'=iA'II(g 2 ; to, to') 



7r lim lim pm(%, x — S) . 
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Knowing that 75 reflects left-handed quark cigcnfunc- 
tions and leaves right-handed eigenfunctions unchanged, 
we find the bare scalar and pseudoscalar correlation func- 
tions in the appropriate limits: 

lim n s (g 2 ) = 2Nfrr 2 lim lim p M (x, x - 5) , 

lim n PS (g 2 ) = -2N f n 2 lim lim p M (x, x - S) . (8) 

The exact expression for the microscopic two point cor- 
relation function, pm(x, y), was found in Ref. [Bj for chi- 
ral ensembles. The connected part is 
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where x = xnp(0)/Y, with p(X) being the continuum 
single eigenvalue spectrum. 

Fixing S as previously explained, the above equa- 
tion (p|) has the limit 



lim pm(x, x — 5) = 



sin[(57rp(0) /Y 2 
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The distinction between the scalar and pseudoscalar 
correlators naturally lies in the sign difference between 
the two. The argument 5 is a small but finite number 
which can be associated with the Minkowski momen- 
tum, g, carried by the quark-antiquark pair. Expanding 
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Eq. ([n]) in small S and inserting this into Eq. (||) , the 
quark correlation functions are found to be 
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valid in the limit of small q. 

Now, we need only perform the sum indicated in 
Fig. 1. The full correlation function, II r (<7 2 ) , is summed 
asEq. ©: 
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Inserting the solutions for each channel, Eqs. (|Tl|), we find 
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2N f C [7T 2 p(0) 2 + £ 2 ] ' 
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for the scalar channel, where C = 3S 4 /[2iY /7 r 4 p(0) 4 ]. For 
the pseudoscalar channel we obtain 



n p V) = 
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q 2 - 2N f C [tt 2 p(0) 2 - £ 2 ] 
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The scalar propagator has a pole at finite q 2 . This result 
is consistent with our expectation of a massive excita- 
tion. (This is simply the scalar mesonic mode, usually 
denoted a.) On the other hand, we expect that pseu- 
doscalar excitations should be identified with the pion 
and, hence, should be massless. Indeed, by virtue of the 
Banks-Casher relation. 



E = 7rp(0) 



(15) 



we see that the pseudoscalar propagator assumes the 
form 
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exhibiting the pole at q 2 = indicative of a Goldstone 
boson. This cancellation is precisely analogous to that 
which occurs when a gap equation, arising in an effective 
field approach (e.g., the Nambu-Jona-Lasinio or instan- 
ton vacuum models), conspires with the g-independent 
part of the pseudoscalar correlator to leave a massless 
mode. Here, yRMT exhibits a corresponding content in 
that a relation between the chiral condensate and the 
single-quark energy spectrum is vital for recovering the 
pion as a massless excitation. Furthermore, the collapse 
of support at zero eigenvalue with the restoration of chiral 
symmetry S mimics the gap equation's loss of validity 
when the condensate vanishes. 

Finally, we note that the scalar propagator of Eq. ( O ) 
implies an effective mass of M 2 = 6£ 2 . Because of the 
ambiguity of scale in our association between eigenvalues 



and momenta this cannot be taken as a precise numeri- 
cal prediction. However, it is clear that the mass scales 
with the strength of the condensate, S, as is common in 
effective field models. 

The present results were obtained with Gaussian dis- 
tributed random matrix elements. However, the demon- 
stration of massless pseudoscalar and massive scalar exci- 
tations relied only on properties of the microscopic spec- 
tral correlators which are known to be universal. Thus, 
we expect that the present qualitative results will apply 
to a greatly extended class of potentials. 

By following the path of an RPA resummation and 
utilizing exact results from random matrix theory in the 
thermodynamic limit, we have demonstrated that evi- 
dence of spontaneously broken chiral symmetry based 
on the Banks-Casher relation is confirmed by the exis- 
tence of Goldstone modes with mesonic quantum num- 
bers. This result indicates that xR-MT inherently con- 
tains phenomena beyond the mere appearance of a chiral 
condensate. It also suggests that the pattern of sponta- 
neous symmetry breaking found in %RMT is the same as 
that occurring in dynamical field theoretic approaches. 
Finally, the underlying connection between the vacuum 
condensate and the single-quark spectrum is further seen 
to be a feature which chiral random matrix theory shares 
with other descriptions of chiral symmetry breaking. 
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